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Abstract 

We consider a general class of electrically charged black holes of Einstein- 
Maxwell-scalar theory that are holographically dual to conformal held 
theories at hnite charge density which break translation invariance ex¬ 
plicitly. We examine the linearised perturbations about the solutions 
that are associated with the thermoelectric DG conductivity. We show 
that there is a decoupled sector at the black hole horizon which must 
solve generalised Stokes equations for a charged huid. By solving these 
equations we can obtain the DG conductivity of the dual held theory. 
For Q-lattices and one-dimensional lattices we solve the huid equations 
to obtain closed form expressions for the DG conductivity in terms of the 
solution at the black hole horizon. We also determine the leading order 
DG conductivity for lattices that can be expanded as a perturbative series 
about translationally invariant solutions. 



1 Introduction 


The holographic correspondence provides a powerful framework for obtaining precise 
results about strongly coupled systems using weakly coupled gravitational descrip¬ 
tions. A key cornerstone is that the dual description of the field theory at hnite tem¬ 
perature is provided by a black hole spacetime. It is a remarkable fact that various 
properties of the thermal system are captured by the properties of the gravitational 
solutions at the black hole horizon. For example, the entropy of the thermal system 
at equilibrium is given by the Bekenstein-Hawking forumla, S = A/AG, where A is 
the area of the black hole event horizon. Another well-known example is provided 
by the shear viscosity, p, which, for certain classes of black hole solutions, is given 
by r; = dvrs, where s is the entropy density (^[^. In this paper, which expands upon 
and generalises the results presented in [^, we explain how the DC thermo-electric 
conductivity can be obtained by solving equations for a non-relativistic fluid on the 
black hole horizon. Moreover, we will see that, unlike rj, our result for the DC con¬ 
ductivity holds in a very general context, being applicable to arbitrary static black 
holes for which the DC conductivity is hnite. The extension to stationary black holes 
will be presented in [^. 

We hrst recall that the DC thermal conductivity, k, is a very natural observable 
to study in holography. Indeed, in the regime of linear response it determines the 
heat current, or momentum how, that is produced after applying a constant external 
thermal gradient. Although this naively appears to be a low-energy quantity it is 
in fact sensitive to the UV physics. For example, if the system is translationally 
invariant, and hence conserves momentum, then k is inhnite. More precisely, there is 
a delta function in the AC conductivity at ca = 0. To obtain a hnite k it is necessary 
to introduce a mechanism for momentum to dissipate. 

For systems with a U{1) symmetry, another natural observable to consider is the 
DC electric conductivity. Momentum dissipation is also required in order to obtain 
a hnite result when the charge density is non-vanishing. More generally, for these 
systems there is a mixing of electric and heat currents and one should consider the 
matrix of thermo-electric conductivities 

j\_/a aT\ f E 
Q )~\aT FT J -(VT)/T 

where J and Q are the total electric and heat current hux densities (which we de- 
hne precisely later) and E and VT are constant applied electric held and thermal 
gradients, respectively. 
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Holographic lattices provide a natural framework^ for studying momentum dissi¬ 
pation. These are black hole solutions whose asymptotic behaviour at the holographic 
boundary corresponds to the addition of spatially dependent sources to the dual held 
theory. Various holographic lattices with sources that depend periodically on just one 
of the non-compact spatial directions have been constructed by solving PDEs in two 
Constructing lattices that depend periodically on additional spa- 


variables 13 “19 


tial dimensions generically requires solving PDEs in more variables, which becomes 
increasingly challenging at the technical level. An important class of exceptions are 
provided by Q-lattices 20,21 which exploit a global symmetry in the bulk leading 
to a problem involving ODEs. Other constructions involving ODEs use massless ‘ax- 
ions’ to obtain sources linear in the spatial coordinates 22 -27 , or use the metric or 


matter helds to obtain helical sources when D > 5 28 31 . A particularly interesting 
application of holographic lattices is that they can lead to novel incoherent metal 


ground states 21,28,32 , insulating ground states 20,21,28 and transitions between 


them 21 


Connections between holographic lattices and superconductivity have 


also been explored in 31,33,34 


For the special case of translationally invariant black holes at zero chemical po¬ 
tential, while the thermal conductivity is inhnite the electric conductivity is hnite 
and can be expressed in closed form in terms of the behaviour of the solution at the 
black hole horizon 35 . For the above holographic lattices involving ODEs, formulae 
for the thermoelectric DC conductivity, also expressed in terms of the black hole so¬ 
lutions at the black hole horizon, were obtained in [^|^[^|^. These results were 
then extended to a class of one-dimensional holographic lattices in 18 , although 


the details were more involved. Given these results it is natural to anticipate that 
similar results can be obtained for all holographic lattices. Here we show that for a 
broad class of holographic lattices, depending on all spatial directions, in general one 
cannot obtain such explicit formulae for the DC conductivity. However, it is possible 
to obtain the DC conductivity after solving a set of fluid equations on the black hole 
horizon. These equations are a generalisation of the forced Stokes equations for a 
charged fluid on the curved black hole horizon, with additional viscous terms arising 
from bulk scalar helds. Recall that the Stokes equations are a time-independent and 
linearised limit of the Navier-Stokes equations for an incompressible huid arising at 
low Reynolds numbers. We will show that the previous results on the DC conduc- 


^ Another approach is to use massive gravity and various interesting results have been obtained 
e.g. [5 -10 . The DC conductivity has also been examined in the context of translationally invariant 
probe branes e.g. 11 12 . In these constructions a finite DC conductivity can arise because the 


delta function is suppressed by 1/N, where N is the number of branes providing the background 
geometry. 
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tivity can all be obtained as special cases where the Stokes eqnations can be solved 
explicitly in closed form. 

The fact that the flnid eqnations which arise are linear and time-independent is 
not too snrprising since we are calcnlating in the regime of linear response and we 
are calcnlating the DC condnctivity. Similarly, the forcing terms are very natnral 
since they arise from the applied sonrces for the electric and heat cnrrents. On the 
other hand only a snbset of the linearised pertnrbation appears in the eqnations on 
the horizon and it is remarkable that the eqnations form a closed system. 

We emphasise that nnlike in the relativistic flnid-gravity correspondence 37 , and 
the associated non-relativistic limit 38,39 , we do not take any hydrodynamic limit in 
obtaining onr flnid eqnations. In the presence of spatially dependent sonrces a natnral 
hydrodynamic limit wonld arise for temperatnres mnch bigger than all other scales, 
inclnding those of the lattice. By contrast onr resnlts are valid for all temperatnres. 
Onr resnlts differ bnt are also reminiscent of the “membrane paradigm' 


more recent worl^which relates solntions of the non-linear Navier-Stokes eqnations 


40 and the 
on 


hypersnrfaces in Minkowski space to obtain black hole solntions (see also 42 -S]). 
We expect that the time-dependent and non-linear versions of onr eqnations will 
play a role in stndying momentnm dissipation for holographic lattices, possibly after 
taking a hydrodynamic limit and this will be reported on elsewhere. Discnssions of 


momentnm dissipation, condnctivities and hydrodynamics can be fonnd in 45 48 


and some recent low-freqnency condnctivity resnlts are presented in 49 


The plan of the rest of the paper is as follows. In section 2 we introdnce the 
holographic model and the class of electrically charged black hole solntions that we 
shall be considering. In section 3 we analyse the linearised pertnrbations, containing 
sonrces for the electric and heat cnrrents, which are associated with the DC con¬ 
dnctivity. We will show how the Stokes eqnations can be obtained by expanding 
the Hamiltonian, momentnm and Ganss law constraints, associated with a radial 
hamiltonian decomposition, at the black hole horizon. The flnid eqnations determine 
electric and heat cnrrents at the horizon in terms of sonrces for the electric and heat 
cnrrents. In tnrn these can be nsed to obtain snitably dehned constant electric and 
heat cnrrent flnxes which are independent of the radial direction and hence give rise 
to the DC condnctivity. If the deformed CFT is living on then the DC condnc¬ 
tivity is a bi(T,d) X bi(T,d) matrix, where bi(T,d) is the first Betti nnmber of We 
emphasise that we provide a precise procednre for calcnlating the DC condnctivity 
of the bonndary theory (the spectral weight of a two-point fnnction) by solving an 


^Ref. also contains a discussion and references to some of the earlier work on fluids and black 
hole horizons. 
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auxiliary set of fluid equations on the black hole horizon. 

In section 4 we analyse some examples. We first generalise our results to an 
arbitrary number of scalar helds which allows us to reconsider Q-lattices. For the Q- 
lattices and also for general one-dimensional lattices, we show that the fluid equations 
can be explicitly solved and we can obtain formulae for the DC conductivity explic¬ 
itly in terms of the black hole solution at the horizon. We also examine holographic 
lattices that can be obtained as perturbative expansions about translationally invari¬ 
ant solutions, including the AdS-RN black brane. We show that the leading order 
DC conductivity can also be fonnd in closed form. We briefly conclude in section 
5 where we pnt some of the main results in a more general setting of general static 
black hole spacetimes. The paper contains fonr appendices, inclnding a discussion of 
the Hamiltonian decomposition of the eqnations of motion with respect to the radial 
coordinate in appendix A. 

2 The background black holes 

We will consider theories in D spacetime dimensions which conple the metric to a 
gauge-held. A, and a single scalar held, 0. The extension of our analysis to additional 
scalar helds is straightforward as we will discuss later. We focus on D > 4. The action 
is given by 

S = J - V{4>) - 4M _ 1 (9.^)2^ . (2,1) 

The eqnations of motion are given by 

( 0 ) ° ’ 

[z(0)Fn = o, 

V20-W(0)-^Z'(0)F2 = o. (2.2) 

The only restrictions that we will make on the fnnctions V (0), Z{(f)) is that V (0) = 
— (D — 1)(D — 2), W(0) = 0 and Z(0) = constant. This ensures that a unit radius 
AdSo solves the equations of motion with 0 = 0 and this is dnal to a CFT with a 
stress tensor, dnal to the metric, a global U{1) current, dnal to A, and an additional 
operator dnal to 0. When Z'{0) = 0, the D-dimensional electrically charged AdS- 
Reissner-Nordstrom black hole also solves the equations of motion and describes the 
CFT at constant charge density. Note that we have set IGttG = 1, as well as setting 
the AdS radius to unity, for convenience. 
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We will focus on a general class of electrically charged static black holes with 
metric and gauge-held given by 

= —UGdti^ + ^ -|- ds^(Srf), 

A = atdt, (2.3) 

where ds‘^{T,d) = gij{r,x)dx^dx^ is a metric on a (d = D — 2)-dimensional manifold, 
Tid, at hxed r. In addition, U = U{r), while G, F, at and 0 are all functions of (r, x*). 
In section we will discuss our main results in the context of a more general class of 
static black hole solutions. 

Asymptotically, as r —>■ oo, the solutions are taken to approach AdSr, with 

U^r'^, F^l, G^G{x), gij{r,x) ^ r^Qijix), 

at{r,x) —)■ /r(x), (j){r,x) —>■ . (2.4) 

The spatial dependence of the boundary metric given by (5(x), ^p(x) corresponds to 
providing a source for the stress tensor in the dual CFT living on M x Similarly, 
/i(x) is a spatially dependent chemical potential for the global abelian symmetry and 
0(x) gives rise to a spatially dependent source for the associated dual operator, which 
we have assumed has scaling dimension A. 

A particularly interesting class of black holes is associated with adding sources 
to CFTs in hat Minkowski space, In this case is topologically M'^. Periodic 
lattices, which have been a focus of study, are obtained by taking the functions G(x), 
gij{x), fi{x) and 0(x) to be periodic functions on If we denote the period in each 
of the spatial directions be Li then for this class of black holes we can, in ehect, take 
Tid to parametrise a torus with periods x* ~ x* -1- Lj. 

The black hole horizon, which has topology is assumed to be located at r = 0. 
By considering the Kruskal coordinate v = t + -f- ... we deduce that the near 

horizon expansions are given by 

U{r) =r(47rT + f/(^)r + ...) , 
at{r, x) = r (x) -f (x) r -|- ... j , 

G(r, x) = (x) -I- G^^^ (x) r -f ... , 

F(r,x) = (x) -I- {x)r + ... , 

(0) I (1) I 

9ij = 9lj +9\jrF... , 

cj) = (j)^^\x) + r(j)^^\x) -1-... , (2.5) 
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with 


G'(o) (2:) = (x) . 


( 2 . 6 ) 


We have added the extra factor of in the leading expression for at for convenience. 
For later use, we observe that electric charge density at the horizon is simply 


PH = ^Z{^)F^^\h = , 

where For the averaged holographic charge density, p, we have 


(2.7) 


P = 


voL 




voL 


d xpH , 


( 2 . 8 ) 


where vol^ = J d'^XyJ—g is the volume of the spatial metric at the AdS bound¬ 
ary. This result follows from the fact that the gauge-equations of motion implies 
dr{y/-^ZF^'^) + di{^/^ZF^'^) = 0 and for the case of non-compact we have as¬ 
sumed any boundary terms vanish. 

This class of black holes includes almost all of the holographic lattices that have 
been constructed to date as special cases. For example, periodic lattices with mod¬ 


ulated chemical potential with non-vanishing zero mode were studied in 13,16,18 


while the case of vanishing zero mode was studied in 15 . Periodic lattices with a 


single real scalar held have been studied in 14,19 . These examples have spatially 


inhomogeneous sources in one direction only. By contrast, the Q-lattice construction 
using two (or more) scalar helds the non-periodic ‘axionic’ lattices 


studied in 22-27,36,51 are homogeneous and the sources can be in any number of 


the spatial directions. Other homogeneous constructions using D = 5 helical lattices 


have been studied in (an additional gauge held is needed to be included to 

cover the examples of l^)- Metric deformations in one spatial dimension were stud¬ 


ied in 17 . Holographic lattices in the presence of magnetic helds have been studied 
in 52 ; the generalisation of our results to include magnetic helds will be discussed 


elsewhere [^. 


3 Perturbing the black holes 

We want to study the holographic, linear response of the black holes after applying 
suitable one-form sources {E, () on for the electric and heat currents, respectively. 
Generalising we incorporate the sources by the addition of terms that are 
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linear in time. Specifically, we consider the following linear pertnrbation: 


5 [ds^) = dg^ydx^dx'" — 2tMQdtdx ^, 

5A = Sa^dx^ — tEidx'^ + tNQdx^ , 

5(t). (3.1) 

Here, 6(j) are all fnnctions of (r,x*), while Ei = Ei{x), Q = We 

demand that E, ( are closed one-forms on Ecj: 

d{Eidx^) = diCidx^) = 0 . (3.2) 


This means that the one forms E, ( can uniquely be written as the sum of a harmonic 
form plus an exact form on In fact it is the harmonic piece that is important 
Later we will see later that the harmonic piece is the important part of the source. 
In the case that is or a torus, for example, we could take an independent basis 
of sources to be the d one-forms Eidx^ (no sum on i) and the d one-forms Qdx^ (no 
sum on i) with constant Ei, Q. 

The functions M,N in (3.1) depend on (r, x*) and we will fix them in terms of 
the background black hole solution via 


M = GU , N = at. 


(3.3) 


This is assumed, in addition to (3.2), in order to solve the time dependence of the 
equations of motion at linear order. In general, this ansatz for the perturbation 
contains some residual gauge symmetry, which, for our purposes, will not need to be 
hxed. 

At the AdSo boundary, as r —>■ oo, we will demand that the fall-off of 5g^y,5a^, 
6(j) is such that the only applied sources are parametrised by {E, (). The asymptotic 
fall-offs of 5g^y,5a^, 5(j) are associated with currents and other expectation values that 
are produced by the sources. The resulting currents will be our primary interest in 
this paper. 

At the black hole horizon, as r —)■ 0, regularity implies that we must have 


Sg„ = U (r) (ift™ (i) + 0(r)) , 6,j„ = t (6,/S’ S) + 0(r )), 

^9ij = ( 3 ^) + 0{r), 6gtr = SgS'’ (x) + 0{r ), 

Sg,. = (I^) - AKi ^ + 0(r), i (SgS (x) + 0(r)) , 

Sa, = (x) + 0(r), Sar = ^ (&'"> (x) + 0(,r)) , 
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with the following constraints on the leading functions of x: 


+ ^g^rr - 2 ^gft = o, ^g^ = ^g 


ti ’ 


ia<“* = iaf. 


(3.5) 


It is worth emphasising that the logarithm terms that appear in (3.4) are a direct 
consequence of the applied sources {E, (). For the scalar held we have 6(j) = (50^°^(x) + 
0{r). 

In the Kruskal coordinates, at leading order in the expansion in r we have 


6{ds^) ~ + 2dvdr{—6g[^"’ + dgl^^) + 25g^^^dvdx'' + Sg^j^dx^dx^ , 

+ 2{vG^°^Ci - 6gll^ + 6gl]'^)drdx^ + ^ (^6g^l'> + 6gll'> - 2 6gi]'^^ dr^ , 

6A ~ da'pdv + — vE^ dx* + — 8a'p)dr . (3.6) 

Note that to obtain the leading order pieces in the perturbed held strength one should 
calculate the held strength hrst and then take the limit r —)• 0. 


3.1 Electric current 

We dehne the bulk electric current density via 

r = ^Z{(8)E^G (3.7) 


When evaluated at the AdSo boundary we obtain J*|oo which is the electric current 
density of the dual held theory as explained in appendix 

In the background geometry we have J* = 0. At linearised order for the perturbed 
black holes we have 


f 


(FG')V2 


GUZ{(8) 



8grt 


d^Qji 


^gjt 

GU 


tMQ 

GU 


djdcLf 


{drSttj + tdrN(j) 




(3.8) 


and we see that the time-dependence drops out because of (3.3). 

The gauge equations of motion given in (2.2) can be written in the form 


di.r = 0, 

drr = dj (y=^Z(0)F^*) , (3.9) 


as well as 


dr {V^Z{<P)E^^) = -d, (^Z(0)F^*) . (3.10) 
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For later use we also note that the perturbation satishes, at linearised order, the 
condition 

(i(4 * 2'(0)F) = 0 , (3.11) 

where k = dt- Indeed this easily follows by writing the components of the d — 1 form 


as 


ik * Z{(I))F =(-l) 


d-2 


id 


... id-ij)F A • • • A dx^'^~^ 


+ 2 (d- 2 )! ^(^i • • • id- 2 jk)\/^Z{(j))F^^dx^^ A ... A dr 


(3.12) 


where e is the alternating symbol with e(l... d) = 1 and then using (3.9). Note 


that in the special case when (^ = 0 we have that fc = is a Killing vector with 


CkF = = 0. It is then very easy to establish (3.11). When 7 ^ 0, fc is no longer 


a Killing vector and furthermore CkF 7 ^ 0. Nevertheless, at linearised order, we still 


have (3.11) as we have just shown. 


3.2 Heat current 

We now dehne the bulk heat current. To do so we want to identify equations of 
motion involving the metric perturbation that have a similar structure to the gauge 
equations of motion. We do this usin^the vector k = dt- The procedure is slightly 
subtle when ( ^ 0 since in this case dt is no longer a Killing vector. We proceed as 
follows. Consider a general vector k which satishes 


= 0, = aV, 


(3.13) 


for some function a. Note, in particular that a Killing vector satishes these conditions 
with a = 0. The conditions ( |3.13[ ) imply that 

= {-R\ + a5l)k^ . (3.14) 

We next write ip = ikA and ikF = dO + 'ip^ with a one-form and 9 a globally dehned 
function. In the special case that CkF = 0 we have dtj: = 0. We now dehn^ the 
two-form G\ 

‘IZ (0) , r,, . 1 




D-2 


-k^rF^^A^ - 


{D-2) 


[{?,-D) d + p>\Z{ct))F^'^. (3.15) 


^Heuristically, one can view this as a Kaluza-Klein reduction on the time direction. Alternatively, 
the analysis is inspired by derivations of the first law of black hole mechanics e.g. 


"^This definition slightly differs from the definition used in 36 . The expression here has the 
advantage that it is globally defined. 
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If we assume that Ck(t> = 0, using the equations of motion (2.2), we can deduce that 


V„G'“' = fa + k- - ^ Z(4,)A,U.(F‘-n ■ (3.16) 

For our setup, with k = dt and working at linearised order, we can choose 6 = —ip 
and we have 

a = - v!‘”(g;iO) - \d. log(GiF)9« 0 . 

(f = —9 = at + 6at , ip = —Ei dx^ + at Ci dx''. (3-17) 


We now define the bulk heat current density via 


Q* = . 


(3.18) 


When evaluated at the AdSo boundary, we show in appendix that Q*|oo is the 
time-independent part of the heat current density of the dual CFT: 


G'-'"v®(Gi“ - 1‘f) = Q'U - tG^''"s/Wdt‘K, 


(3.19) 


Here t**, j* are the expectation values of the holographic stress tensor and current 
vector, with e.g. J*|oo = The precise combination that appears on the 


left hand side in (3.19) is the operator that is sourced by —tQ and is, by dehnition. 


what we call the heat current density. Notice that in the case when the holographic 
lattice has no spatially dependent sources for the metric this reduces to the standard 
expression t** — /ijh The time dependent piece on the right hand side is associated 
with the static susceptibility for the heat current two point function (see appendix C 

of 1^). 

In the background geometry we have Q* = 0. At linearised order for the perturbed 
black holes we have 




Q3/2JJ2 

FV2 


^^91 d, 


ddjt 

GU 


dj 


dg 


rt 


GU 


atT 


(3.20) 


W = 0, 

drQ^ = d, , (3.21) 

as well as 

a, (y^SG") + a, (y= 5 G‘‘) = ffa + - ^z(,),)F 6 >d . (3.22) 
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Note that (3.22) includes an equation for the background as well as the linearised 
perturbation. We also record here that 






-|- difdciij —k I 

(3.23) 


and we note that the time dependence drops out because of the conditions (3.2) and 


(3.3). 


Finally, following the discussion for the electric currents, with k = dt we conclude 

d{ik * G) = 0 . 


that (3.21) implies 


(3.24) 


3.3 Currents at the horizon 


We now obtain expressions for the electric and the heat current densities by expanding 
at the black hole horizon. We hnd: 




j;„) = J'|„ = ' + Ej) - 

'3(0) = '3‘Ih = . 


From the hrst equations in (3.9) and (3.21) we immediately obtain 

diJlo)=0, diQ\Q-^ = 0, 


(3.25) 


(3.26) 


which give two equations for a subset of the perturbations at the horizon. We can 
obtain a closed system of equations, which are the generalised Stokes equations, by 


considering the second equation of (3.21). We explain how this can be achieved in 


appendix The same system of equations can also be obtained, in a more illuminat¬ 
ing manner, by evaluating the Hamiltonian, momentum and Gauss law constraints 
on the black hole horizon, as we now discuss. 


3.4 Constraints at the horizon 

We carry out a Hamiltonian decomposition of the equations of motion using a ra¬ 
dial decomposition in appendix The momentum constraints and the Gauss-law 
constraints, can be written in the form H'' = C = 0 where 

((_/i)-V2^/..) + , 

G , (3.27) 
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with 


- K , 
vr^ = V^ZF^^Pup , 

t^4> = —'J—hrFdycj). (3.28) 

Here n is the unit norm normal vector, h^i, = — riprii, is the induced metric, 

Kpu = \^nhpu is the extrinsic curvature and K = In addition Dp is the 

Levi-Civita connection with respect to hp^, bp = hp'^At, and fpu = dphy — d^hp. 

We want to analyse these constraints for the perturbed metric on a surface of 
constant r, near the horizon, and then take the limit r —)■ 0. In local coordinates we 
have n = Ndr where N is the lapse function. We immediately notice that 

, (3.29) 


and hence the Gauss-Law constraint G = 0 is simply dpJ^ = 0. Evaluated at the 
horizon we obtain = 0 as in (3.26). Turning to the momentum constraint, which 
we discuss further in appendix]^ we hnd that evaluating Ht = 0 as an expansion at 
the horizon gives = 0 as in (3.26). Next, Hi = 0 evaluated at the horizon gives 

the extra equation mentioned above which, combined with = 9jQ(g^ = 0 gives 

the Stokes system of equations which we summarise in the next subsection. Finally, 
we note that the leading order term of the Hamiltonian constraint, which is explicitly 
given in (A.9), also gives the condition 9iQ(g^ = 0. 


3.5 Generalised Stokes equations at the horizon 


We can now summarise the closed system of equations that we have shown a subset 
of the linearised perturbations must satisfy at the black hole horizon. The black hole 
horizon is as in (2.3), (2.5) and (2.6). The perturbation at the horizon is given as in 
(3.4),(3.5) and it is illuminating to now introduce the following notation: 




w = Saf'\ 


p = -AtiT lnG(°). 


(3.30) 


These d -|- 2 unknowns satisfy the following d -|- 2 linear system of partial differential 
equations: 


ViF = 0, (3.31) 

Vi{Z^^'>V^w) + F Vi = -Vi(Z(°)F*), (3.32) 

-2 VWpUj) - Z^^^af^VjW + = AnT Q + , (3.33) 
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where the covariant derivatives V in this snbsection are with respect to the metric, 
on the black hole horizon and all indices are being raised and lowered with 
this metric. The hrst two eqnations are simply = diJio) = 0, where 

are the heat cnrrent and electric cnrrent densities at the horizon, respectively: 

0 ( 0 ) = , 

4) = v^Sw)Z'°’ (%“ + + Ej) . (3.34) 

It is helpfnl to note that in the third eqnation we can also write 

2 = V\- + RjiV^. (3.35) 


We emphasise that by evalnating the constraints at the horizon we obtain a system 
of eqnations for a subset of the linear pertnrbation, namely, ^5'^^ ^ind we 

obtain a closed system for this set. Fnrthermore, the eqnations we have obtained are 
a generalisation of the forced Stokes eqnations for a charged fluid on the curved black 
hole horizon. Indeed in the special case of electrically neutral black hole horizons 
with = w = E = 0 and in addition constant the equations are simply the 
Stokes equations with fluid velocity Vi, pressure p and forcing term given by the closed 
one-form 4:7tT(. The curvature of the horizon gives rise to an extra viscosity term 


as in (3.35). In the general case we have a charged fluid with scalar potential w and 


an additional forcing term given by the closed one-form E. It is also interesting to 
note that the scalar field is giving viscosity terms of the form We 

will see in section 4 how these extra terms play a direct role in determining the DC 
conductivity. We emphasise that we have not taken any hydrodynamical limit in 
obtaining these equations. 

We now establish a number of interesting properties of this set of equations. 


Firstly, by taking the divergence of (3.33) and using (3.31), (3.32) we obtain the 
“pressure Poisson equation” 

VV = Vj + E^) + dvrTC^' - j • (3-36) 

For a compact horizon and given background data, the pressure term is uniquely 
specihed by Uj, w, Ej, Q. 

Second, we multiply (3.33) by Vj from the left, and then integrate over the horizon. 


and use (3.31),(3.32) to obtain 


^ {Vw + Ef + Vi Vj\ 

= [ d^x (Q(o)O + Jlo)E^) (3.37) 
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In the case of non-compact horizons, we have assumed that possible boundary terms 
vanish. Observe that the left hand side is a manifestly positive quantity and this is 
related to the positivity of the thermoelectric conductivities, which we discuss later. 


Third, we consider the issue of uniqueness of the equations (3.31)-(3.33). If we 


have two solutions then the difference of the functions, which we again write as 
{vi, w,p), will satisfy the same equations but with vanishing forcing terms, Q = Ej = 


0. From (3.37) we immediately conclude that 


= 0 , = 0 , = 0 . 


(3.38) 


We also have = 0 from (|3.32|) and Vjp = 0 from (|3.33|). We conclude that the 


solution space of equation (3.33) is unique up to Killing vectors of the horizon metric. 


with p, tc constant. We then have 5grt = (dvrT)”^ plus a constant. This result 

agrees with the intuition that one should be able to boost along the orbits of Killing 
vectors to obtain a solution with momentum at the horizon. 

Fourth, we observe that when (i7, Q are exact forms, (i7, C) = (de, dz) with e, 2 ; 
globally dehned functions on we can solve the equations ( |3.31| )-( |373^ by taking 
w = —e and p = AtiTz, plus possible constants, and u* = 0. We observe that this 


solution gives no contribution to the current densities (3.34) at the horizon. We will 


see later that this solution gives no contribution to suitably averaged currents at the 
AdS boundary and hence no contribution to the DC thermoelectric conductivity i. e. 
the DC conductivity is determined by the harmonic part of E and C,. A basis for 
the non-trivial part of these sources is thus given by a basis for the hrst cohomology 
group of Erf. 

Fifth, we point out that the fluid equations that we have obtained can be obtained 
by varying the following functional: 


L= / ^ 


+p{ViE) + -f Vw)' 


+ AttTCiE + Z^^'^Ei + V*w) + -Z^^^EiE^ 


(3.39) 


and we remind the reader that the covariant derivative V in this subsection is with 
respect to the metric Varying with respect to the pressure, which is a Lagrange 
multiplier, gives the incompressibility condition. Varying with respect v and w then 
gives the remaining Stokes equations. It is also interesting to note that if we vary with 
respect to Ei and Q then we get the currents at the horizon and Q(o)’ respectively. 
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On shell we therefore can deduce, for example, that 


51, 


( 0 ) 


5Q^ 


( 0 ) 


50 


5E, 


(3.40) 


This is a kind of Onsager reciprocal relation for the currents at the horizon. After 
considering the current fluxes, to be described in the next subsection, we obtain 
Onsager relations for the DC conductivitities. 

Finally, we comment on the fact that, locally, the sources can be eliminated from 


the Stokes equations (3.31 )-(3.33). Indeed since the sources E,C. are closed, locally 


we can write E = de, C, = dz and after dehning w = w + e, p = p — AnTz we have 


= 0, 

+ O Vi = 0, 

-2 = 0 . 


(3.41) 

(3.42) 

(3.43) 

It is important to emphasise that now w and p are not globally dehned functions on 
the black hole horizon. For example, if the horizon was a torus with x* = a:* + Li, and 
the source E = cdx^ for some constant c, then w would satisfy the twisted boundary 
condition w{x^ + Li) = w + cLi. Note that this would give extra contributions to 


(3.37). It is therefore most natural to work with the formulation with sources, and 


we will do so in the sequel. 

It is worth noting, however, that the sources can also be removed, locally, from the 
full linearised perturbation. Indeed, suppose we carry out the gauge transformation 
A —)■ A = d{te) + B and in addition change the time coordinate via t = t(l — z) + C, 
where B, C are functions independent of the time coordainte. We then hnd that at 
linearised order we obtain the same perturbed ansatz with vanishing sources and 

+ 2f/ Gz, 5gi^ = 6gu - UGdrC dgn = 5gu - UGdiG 
5ai = Sat + e — atz, Sar = Sur + drB, Sat = Sat + diB (3.44) 

We choose B, G to vanish suitably fast at the AdS boundary and at the horizon we 
choose i? = lnr/(47rT)e + ... and C* = — lnr/(47rT)z + .... Evaluating at the horizon 
we see that in the new coordinates and gauge we have induced w ^ w, p ^ p and 


Sgi? Sgtt 


( 0 ) 


+ 2G^°h. 


3.6 The DC thermoelectric conductivity 


For a given set of sources {E,() we can solve the Stokes equations (3.31)-(3.33) 
at the black hole horizon and hence obtain expressions for the electric and heat 
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current densities at the black hole horizon. From this data we would like to deduce 
something about the current densities at the holographic boundary as a function of 


(EX)- The radial dependence of the current densities are given by (3.9) and (3.21). 


For some simple special classes of black hole solutions the current densities 
are independent of the radial coordinate. This occurs for the Q-lattices and the 
holographic lattices that depend on just one of the spatial dimensions, for example. 
However, for general classes of black holes J*, Q* will depend on the radial direction. 

On the other hand, remarkably, we can always dehne “current flux densities” J*, Q® 
which are independent of r and hence we can obtain the associated DC conductivity. 
To consider a concrete example, we assume that we are in D = 4 with a periodic 
holographic lattice and or a two-torus. In particular, the lattice deformations 


G{x), gij{x), /i(x) and (j){x) in (2.4) are all periodic functions of the spatial coordinates 
X® with period Lj. We dehne the following current hux densities 




J^ = -r J^dx^ . 
Li J 


(3.45) 


where X and is the current hux density through the X 2 and Xi planes, respectively. 


We dehne Q* in a similar way. We can then immediately deduce from (3.9), (3.21) 


that drJ'' = drQ^ = 0, which is simply Stokes theorem in the bulk. Notice that 
J®, Q® are also independent of the x® coordinates and hence they are just constants. 
Similarly in D = 5 with a periodic lattice on we would dehne the following 
constant current hux densities 


Ji = 


1 


L 2 L 3 


dx^dx^J^ 


J2 = 


L 3 L 1 


dx^dx^j"^ 


= 


L 1 L 2 


dx^dx^J^ , 
(3.46) 


The current hux densities of the boundary theory are, by dehnition, J®, Q® evalu¬ 
ated at the AdS boundary. In order to evaluate the DC conductivity matrix we want 
to relate these to constant sources Ei, Q at the AdS boundary via 



a 






Td®^' Th®^' 



(3.47) 


We have just shown that J®, Q® are constant and so their value at the AdS boundary is 
the same as at the black hole horizon, and that in turn these are hxed by the closed 


forms EX evaluated at the horizon by solving the Stokes equations (3.31 )-(3.33). 
Continuing with the case that or a torus, we can take an independent basis 

of sources to be the d one-forms EXx^ (no sum on i) and the d one-forms C^dx® (no 


sum on i) with constant EiXi and this dehnes the DC conductivity matrix (3.47). 
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Note that an equivalent way to characterise the constant sources Ei, Q at the 
AdS boundary is to write a general closed form source on or a torus, as 

E = Eidx^ + de, where e is a periodic function and then extract Ei by integrating 
E over an appropriate basis of one-cycles. It is worth emphasising that in the para¬ 


graph following (3.38) we noted that it is only the harmonic part of the sources that 


contribute to the currents at the horizon and hence this procedure gives the same 
DC conductivity. 


After substituting (3.47) in (3.37), we can now explicitly see that the positivity of 


the left hand side of equation (3.37) implies that the thermoelectric DC conductivity 


is a positive definite matrix. Continuing on from the discussion following (3.39) we 
can deduce that the thermoelectric matrix is symmetric. 

We will discuss the DC conductivity matrix when is not or a d-dimensional 
torus in the next subsection. 

We have presented a procedure for calculating the DC conductivity of the bound¬ 
ary theory in terms of a calculation at the black hole horizon. One might wonder if 
the calculation could also be done at any constant radial hypersurface. In fact this 
cannot be done since evaluating the constraint equations on a constant r hypersurface 
with r tA 0 will not lead to a closed system of equations for a subset of the linear 
perturbation and hence we cannot obtain the current fluxes. 


3.7 Perspective using forms 


We have focussed on a particular class of black holes given in (2.3), with a single black 


hole horizon, and with perturbation given in (3.1) with (3.2), (3.3). We have also fo¬ 


cussed on black holes for which Yd is topologically either or a d-dimensional torus. 
In this section we briefly discuss the DC conductivity calculation using the language 
of forms, which illuminates some global issues as well as revealing generalisations for 
Yd with other topologies. 

One key point is that the two-forms E and G for the perturbed metric satisfy, at 


linearised order, the following closure conditions (see (3.11), (3.24)): 


d{ik* Z{(j))E) = 0 , d{ik*G)=0, 


(3.48) 


where k = dt- These conditions are valid without k being a Killing vector, but instead 


satisfying the weaker conditions in (3.13). 


We now observe that for any two-forms satisfying (3.48) we can define a natural 


set of current fluxes. Specihcally, at the deformed AdSo boundary and at hxed t, we 
let Ca, a = 1,..., bd-i{Yd), be a basis of d — 1 closed cycles on Yd, where bd-i{Yd) is 
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the {d — 1) Betti number of We can then dehne the current fluxes through these 
cycles via 

r = - f tk*Z{^)F, f tk*G. (3.49) 

^ Ca Ca 

We now consider a d-dimensional surface S in the bulk spacetime which has boundary 
Ca at the AdS boundary and possibly another boundary at a black hole horizon. Then 


since the integrand in (3.49) is closed, we deduce that these current fluxes are also 


equal to their values at the black hole horizon. If the cycle Ca is contractible in the 
bulk, then the current flux would necessarily have to be zero. 

For the special cases of periodic lattices in D = 4 and 5 spacetime dimensions, 
for which is topologically and respectively, using (3.12) we immediately 


see that the dehnition (3.49) agrees with the dehnition^ given in (3.45), (3.46) after 
choosing an obvious basis of one and two cycles, respectively. For these cases, the 
number of current fluxes is the same as the dimension of However, this is not 
the case for more general In H = 4, for example, we can envisage black holes in 
which S 2 is a Riemann surface with genus S' > 1, and it is possible to dehne 2g current 
huxes. There are many more possibilities for solutions in H = 5. We also note that 
when Erf is a sphere, which is relevant for solutions associated with deformations of 
global AdS, these current huxes are all trivial since 6^-1 (5''^) = 0. 


The above comments were based on general two-forms satisfying (3.48). A sec¬ 


ond key point in our derivation of the DC conductivity is that the two-forms were 
constructed with specihc source terms parametrised by the one-forms E, (. For the 
class of solutions that we considered we assumed there was a single black hole hori¬ 
zon with the same topology as the spatial boundary of the deformed AdS space. 


In order to satisfy (3.48) it was necessary to take the one forms EX to be closed 


one-forms on and independent of the radial coordinate. Corresponding to the 
basis of {d — l)-cycles, Ca, we can dehne a basis of harmonic one-forms, on 
by Poincare duality. We can then write E = EaX + de and ( = Ca0“ + dz with 


constant Ea, Ca- Recalling the discussion in the paragraph following (3.38), in solving 


the Stokes equations at the horizon in order to obtain the currents at the horizon 
only the harmonic part of the sources, EaX, CaX important. We therefore can 
relate the hd-iiXd) independent constant source terms to the bd-i(T,d) current huxes 
at the horizon, after solving the Stokes equations, and hence to the bd-iCZd) current 
huxes at the AdSr, boundary. This procedure give rise to thermoelectric conductivity 
matrices (t“^, a°'^ and all of which are 5i(Sd) x 61 (S^) matrices, where we 

used the fact that bd-iiXd) = 61 (E^). 


Note that in (3.451, (|3.46[) we have divided by suitable Li in order to obtain current flux densities. 
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One can ask if this aspect of the formalism can be adopted to more general classes 
of black holes in which there are multiple black hole horizons (an example of such 
a solution, but without spatially dependent sources, is given in [^). This would 
require identifying suitable source terms ( that depend on the radial direction 
while still maintaining the condition (3.48). We return to this point in section]^ 


Finally, recalling (3.37) we note that we can write 


d'^X {Q\o)Ci + '^( 0 )-^*) 


(4 * Z{(j))F) A 4 + (4 * O') A C , 


= - / {ik* Z{(t))F) - {ik*G), (3.50) 

Jc(^) Jc(0 

where in the hrst line we are integrating over any surface at constant r and t. In the 
second line 4^®^ and are (d — 1) cycles, unique up to homology, that are Poincare 
dual to the closed one-forms E and C,. By dehnition the right hand side is thus the 
sum of the current fluxes through the cycles and repsectively. 


The positivity of the left hand-side, which we obtain from (3.37), is associated with 
the positivity of the thermoelectric DC conductivity. 


4 Examples 

In this section we examine some special examples of holographic lattices for which 
we can solve the fluid equations on the horizon and hence obtain expressions for 
the DC conductivity in terms of the behaviour of the black hole solutions at the 
horizon. We hrst discuss how extra scalar helds manifest themselves in extra terms 
in the huid equations at the horizon and then use this to study a general class of 
Q-lattices. We next analyse general holographic lattices that depend on just one 
spatial dimension. Finally we examine holographic lattices that can be obtained as 
a perturbative expansion about the AdS-RN black brane. 


4.1 Extra scalars and Q-lattices 


For simplicity we derived the Stokes equations for the model given in (3.33) which 
involved a single scalar held. However, the generalisation to extra scalar helds, which 
can parametrise a non-trivial target space manifold, is straightforward. Specihcally, 
if we replace (2.1) with several scalars, 0^, with the functions V, Z depending on all 
of the scalars and the kinetic-energy terms generalised via 

-la.#-" ^ -^euws-p'd ^-', (4.1) 
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then this leads to the Stokes equations as before, with the only change in (3.33) given 
by 


With this result in hand we can now obtain previous results for the DC conduc¬ 
tivities for Q-lattices 21,36 . They key feature of the Q-lattice is that it exploits a 
global symmetry in the bulk to construct the black hole solutions. In the present 
context we assume that the model admits n global shift symmetries of the scalars: 


-)■ 


jjia ^la 


(4.3) 


with a = l,...,n. For example, if we had a single complex scalar held with a 
global U{1) symmetry, then the associated shift symmetry of this type is obtained by 
parametrising the scalar manifold locally with the modulus and phase of the complex 
scalar held. This gives rise to a periodic lattice. Another example, is a massless ‘axion’ 
held with only derivative couplings. Note that the function Z must be independent 
of these 

The spatial coordinates are taken to parametrise or possibly a torus. The 
black hole solutions are then constructed based on an ansatz in which the scalars 
associated with these shift symmetries take the form 

^ qIc. ^ 


everywhere in bulk with C a constant n by d matrix. For simplicity of presentation we 
assume that all spatial coordinates are involved and hence the DC conductivity in all 
spatial directions is hnite. The metric, the gauge-held and the remaining scalar helds 
will depend on the radial direction but will be independent of the spatial coordinates 
X*. The metric on the black hole horizon is hat and in addition, and a® 

are all constant. 


After these remarks, the huid equations (|3.31|)-(|3.33|) are solved with v\ p and w 

(4.5) 


all constant on the horizon. The huid velocity is given by 


= attt {v~^y 


P 


0 + 


with constant Ei, Q and we have dehned the dx d matrix: 

Vi. = Gi j 

L Irv-i J-an t j 


(4.6) 


Furthermore, the averaged charge density, p, dehned in (2.8), and the entropy density, 
s, are given by 


p = ph = 




s = ■ 


(4.7) 
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The current densities are independent of the radius and are given by their 

horizon values: 


T = 


47r 




Q‘=4nTs(V-'y’ + . 


(4.8) 


The DC conductivities are thus given by 


a 




a 






s j,- 47rp^ 

= ^^( 0 ) + — 

= 47rp {V-^Y 

=47rTs {V-^^ . 


{v-y , 


1 


(4.9) 


Note that the conductivity when Q = 0, ctq=o = (t — Tan is given by 

s ^( 0 ) .. 

= ( 4 - 10 ) 

A hnal point worth emphasising is that the origin of the matrix V appearing in 
the hnal expressions arises from the extra terms involving the scalars in the Stokes 
equations, underscoring the signihcance of the latter. 


4.2 One-dimensional lattices 

We now consider a class of black hole solutions with metrics on the horizon that break 
translations in just one of the spatial directions. As special sub-cases we will recover 
the results for the inhomogeneous lattices with varying chemical potential studied 

Recently formulae for the DC 


m 


18 as well as the helical lattices studied in 29 


conductivity for a scalar lattice were obtained in 19 in terms the behaviour of the 


solution at the black hole horizon as well as sub-leading terms. We improve upon 
those results by providing a new formula that depends just on the solution at the 
horizon. 

We assume that the horizon geometry depends on the spatial coordinate x and 
is independent of the remaining d—1 spatial coordinates which, for dehniteness and 
without loss of generality, we take to parametrise a torus. The moduli of this torus 
can depend on x. For simplicity we restrict our considerations to metrics on the black 
hole horizon of the form 


ds\ = qY dx^dx^ = jIx) dx^ -|- ds]^_i (x) , 


(4.11) 
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where ds‘^_i = gabdx°'dx^ is a flat metric on the torus. We now solve the relevant 


system of equations (3.31)-(3.33). The incompressibility condition (3.31) is solved by 


taking the non-vanishing components of v* to be 


(4.12) 


with Qd-i the determinant of the d — 1 dimensional metric on the torus and vq a 
constant. The non-trivial component of the current density is which must be a 
constant, and we have 


7 


1/2 




'^{0) ~ H-- Vq -\- Ex 


1/2 

9d-i 


(4.13) 


With a little effort we can now write the Stokes equation (|3.33|) in the form 
2 Vo dx dxQd-i^'j -Y Vo 


+ - dxP = -dvrT C. , 


1/2 

9d-i 


(4.14) 


where we have defined 
1 


^ = b?—^TT 72 [i^x'^rigd-iY + dxgabdxgcdg’"" 
2(75'd-i) / 


h9d-i)Y^ 




1/2 

9d-i 


(4.15) 


with gab{x) the metric components for Equations (4.13) and (4.14) may now 

be used to fix the functions w and p. Since these are periodic functions, we must 
have that the expressions for dxW and dxP have no zero modes on the torus and this 
imposes constraints on and Vo- A simple way to establish these constraints is 
take an average of the two equations. In fact doing this completely Axes J(q) and vo 
in terms of the sources. Indeed, if we define 






(4.16) 


where Li is the period of the lattice, we obtain 


Jm = ^ I 


y+ tore. 


9l\ 


no = ^ I 47rT Cx 


7 


1/2 




+ Ex 


9T\ 


(4.17) 
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where 





& 


(4.18) 


The expression for the heat current at the horizon is simply = 47rT uq- Now 
for these one-dimensional lattices the electric current and heat current densities are 
independent of the radial direction and so we have deduced their values at the AdS 
boundary. Thus we can immediately extract the thermoelectric conductivities in the 
X direction and we hnd 


1 /■ _ _ 47r 

a = — Y, a = a = — 


X 


X 


1/2 ' 
9d-i 


K = 


X 


7 


1/2 


1/2 

9d-i 




(4.19) 


Observe that the hnal result for the conductivity is invariant under reparametrisations 
of the X coordinate, as it should be. 

We can also write the electrical conductivity in the form 

2 / , \ -1 


— O'Q=0 



7 


1/2 


d-A'”/ 


(4.20) 


where <Jq=o = a — Tan is the conductivity when Q = 0 (as opposed to C = 0) 
and is given by 


<^Q=o 


7 


1/2 


-1 




(4.21) 


Notice that the second term in (4.20) vanishes if = 0. 

We hnish by indicating how to obtain some previous results. For the one-dimensional 
lattice of H = 4 Einstein-Maxwell theory given in 18 we simply need to use the 
translation given by 


7 = S(a;) dsl = S(x) = H^\x ), 


(4.22) 


in order to obtain the expressions for a, a, a, k given in 18 . Similarly for the helical 

(4.23) 


lattice of pure D = 5 gravity studied in we should use 

Y = hi, dsl = rl (e^“+ ul + e"^“+ ul) , 
where two of the three left-invariant one-forms for Bianchi VIIq are given by 

U 2 = cos kx dy — sin kx dz, = sin kx dy + cos kx dz . (4.24) 
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It is straightforward to show that V = sinh^ 2a+/(r^h+) and hence recover the 


formula for k in the x direction given in 29 . Finally charged helical lattices in D=5 


gravity with two gauge-helds coupled to a scalar were studied in 30 . Setting the 
second gauge-held to zero, we can compare our results by setting 


7 — Cl, ds2 — C2UJ2 + C3CJ3 


A short calculation shows that in the x direction we have 

C 2 C 3 47rp2 


a = 


L/q 03 Z/Q 


^1/2 ■ P(C2-C3)2 S 

C2C3 


+ 


a = a = 


-Attp, 


C2C3 


K = 


-AtisT , 


(4.25) 


(4.26) 


fc2(C'2-C3)2 fc2(C'2-C3)2 

where s = 47 r(CiC 2 C 3 )^C and p = {€ 1020 ^)^/“^(see (|2.8[)). The expression 


for a agrees with (SOl and the expressions for a, a and k are new. 


4.3 Perturbative lattices 

We now consider the case of a periodic lattice that is constructed as a perturbative 
expansion about the electrically charged AdS-RN black brane solution with a hat 
horizon. As we have noted before since everything is periodic in the spatial directions, 
in ehect, we can take to be a torus. If A is the perturbative parameter, then at 
the black hole horizon we will assume that we can write 

g(o)ij = g ^ij + + ■ • •, c*'^^ = /(o) -I- A /(i) + • • •, 

= a -|- A 0(1) = '0(o) + '0(i) + ''' ) 

= ^( 0 ) + ^( 1 ) + ■ ■ ■ ) ( 4 . 27 ) 

with a, 'ip(o), /(o) and g being constant and the sub-leading terms are periodic 
functions of, generically, all of the spatial coordinates xh Note that the entropy 
density and the electric current density on the horizon are given by 

5 = 4 ^/^ PH = ag'^l\ ( 4 . 28 ) 

For the Ricci tensor and Christohel symbols we have the expansions 

-R(o)ii = Rfj^ + • • •, 

r}» = Ar 7 ‘ + A^rf ‘ + ■ ■ ■. (4.29) 
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It turns out that we can solve equations (3.31 )-(3.33) perturbatively in A using 
the following expansion: 


^(0) + y + ^h) + • • • ’ 


A2 

1 


^ ^ A -’ 


P = ^P(l) +P(2) + ••• . 


Expanding (3.31)-(3.33) in A, we hnd at leading order that 

di V(o) = 0 , 


□ — 0 , 


(4.30) 


(4.31) 


where □ = 6'^Widj and we deduce that are just constant on the torus. We proved 
earlier that for the full non-linear problem, in the absence of horizon Killing vectors, 
there is a unique solution to the Stokes equations. Therefore, it must be the case that 
these integration constants are hxed at higher orders in the perturbative expansion, 
and this will be conhrmed shortly. 

At next order in the expansion we hnd 

9iVa) + ^ vLs = 0 , 


9 ^^(0) 


^(0) 

-1- n(o) di = 0 , 


njo) + ad^w^i) - a>(i) = 0 , 


where 






(r('^)L = ^ {dsh^^^'k + dkh^^'>\ - 


(4.32) 


(4.33) 


and all indices are raised and lowered with 6. By considering the pressure Poisson 
equation (3.36) and using VjR^i = h Vji?, we deduce that 


-aw^i)+p^i) = (□ , 


(4.34) 


witht(;(i) = (□ ^(0) from ( |4.32[ ). Combining these results we can obtain 

an expression for in terms of ^(0)- 


^(1) = ^(i)i ^(0) ’ 


(4.35) 


with 


= -D-i - a* . (4.36) 

The function is dehned up to a constant on a torus. The associated constant 

for will be hxed at third order in the perturbative expansion and it does not ahect 
the leading order DC result. 
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We now integrate equation (3.33) to find 

J j 9('oj“ 3)™ + J 9('oj‘9jP+ J 9%^ 

= {J slifj + y g;/"z'“>ar>) , (4,37) 

where we have taken Ei, Q to be constants and J is dehned to be the average over a 
period: 

r 1 


■H- 


Z/l . . . Ln 


dx^... dx^ . 


(4.38) 


We next expand equation (4.37) with respect to A and keep the A° pieces. Using 


terms of v)n\- Indeed we deduce that 


(4.34), (4.35), after some work we hnd that the left hand side can be expressed in 

(4.39) 


40) 


A ‘^LjiVQ = 47rT Q + a U,-, 


where L is a matrix that only depends on the background data given by 

=\^g-^ j d^N\ + i 

+ •^^4^(0) / ®(i) + A^ / di^!(fi)djip , 


(4.40) 


with N as given by (4.36). Notice, in particular, that the integration constants 


associated with □ ^ drop out since the relevant terms are all covered by an extra 
spatial derivative. Thus at leading order we have 


{L-^r {WCj + aE,) , 


(4.41) 


and 


r\H-pHv\ Q%^Tsv\ (4.42) 

Recalling the dehnition of the radially independent current flux densities given in 


(3.45),(3.46), we hnally obtain the holographic current flux densities in terms of E, 

X ^ pv ^, Q* ^ TsE , (4.43) 


where we used (2.8). Thus we can determine the leading order behaviour of the 
conductivities: 

-i4vrp2 


K = L ^AtisT , a = a = L , a = L 


-1, 


(4.44) 
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We observe that K{aT)~^ = s^/p^, which corresponds to a kind of Wiedemann- 
Franz law. In addition we note that the thermal conductivity at zero current flow, 
K = R — Taa~^a, appears at a higher order in the expansion: k ~ A°. It is interesting 
to compare these results to the discussion in [^. We similarly hnd that aQ=o also 
appears at a higher order in the expansion. 


4.3.1 Perturbative one-dimensional lattices 


We conclude this subsection by discussing the special case of one-dimensional pertur¬ 
bative lattices and hence make contact with the results of section 14.21 We hrst notice 
that for an arbitrary periodic function of a single coordinate x we have 


□ ^dxF = c+ / dxF{x)—x / F, 


(4.45) 


where c is an arbitrary constant and in the last term J refers to the average integral 


over a period in the x direction, as given in (4.16). Hence 


J Fdx{n-^dxF) = J F 


(4.46) 


Recalling (4.11) we next assume that 


dx^dx^ = 6'j dx‘^ + 5gab dx°'dx ^, (4.47) 

where dj and 6gab only depend on x. The only non-zero matrix element of the matrix 
L is given by 


A Lxx g 


-2 


^ dxdgab dxdg^^'" + ^ dxdg^a dxdg\ ) / {dx'ipf 


+ 




0 ( 1 ) 


where indices are raised with In obtaining the above result we used that 

Rxx = -^dldg^a, Rab =-^dldgab- 


In the notation of section 14.21 we have 

7 = -h A (57, gab = g dab + ^ dgab , 


(4.48) 


(4.49) 


(4.50) 


and the rest of the functions are expanded exactly as in (|4.27). After this identihcation 

A = —A— Q Lxx . (4.51) 


we hnd that X as dehned in (4.18) takes the form 


J 9 Rxx • 
9 Ro) 
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It is then straightforward to see that the leading order expansion in the lattice 


strength of the DC conductivities give in (4.19) agree with the perturbative results 


given in (4.44) when restricted to lattices that depend on one spatial coordinate only. 


5 Discussion 

The main results that we have obtained in previous sections apply in a more general 
setting as we now discuss. Specihcally, we consider the following ansatz for a general 
class of static solutions 


= gudt^ + A = atdt, 


(5.1) 


where gu^atifl) and the metric ds‘^{MD-i) are all independent of time and are just 
functions of the coordinates on 

The spacetime may have various types of asymptotic boundaries, but our primary 
interest is when there is an AdS boundarjj^ In this case we can introduce a local 


radial coordinate and then impose the same boundary conditions as in (2.3),(2.4), 


corresponding to the CFT living on deformed with various spatially dependent 
sources. The spacetime may have one or possibly more black hole horizons (examples 
have been discussed in j^). Near each black hole horizon we can again introduce 
a local radial coordinate and then demand that the metric has the behaviour that 


we gave in (2.5). Note that we do not assume that the topology of the black hole 


horizons are all the same, nor do we assume that that have the same topology as 
We now consider the following linear perturbation 


6 {ds'^) = Sg^ydx^dx'" + 2tguCadtdx'^ , 
SA = 5a^dx^ — tEadx^" + tatCadx°‘ , 

50 , 


(5.2) 


where the one-forms E, ( are now dehned on Md_i (not just on as before) and 
are still taken to be closed. In addition 5(7^^, 6a^ and 50 are all independent of t. It 
is an interesting fact that at linearised order in the perturbation (and independent of 
any boundary conditions) we still have the key results 


dik * Z{(l))E = 0, dik * G = 0 , 


(5.3) 


where k = dt- Note that k still satishes the conditions (3.13), and we have dehned G 


as before in (3.15) 


®We can also consider other holographic boundary conditions, for example, asymptotically Lif- 
shitz, or even asymptotically flat boundary conditions. 
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We now consider the boundary conditions on the perturbation. At the holographic 
boundary E,( approach closed one-forms ^ind we impose that these 

are the only sources deforming the CFT. Similarly, for each black hole horizon the 
perturbation behaves in local coordinates as in (3.4), with E,( again approaching 
closed one-forms on each horizon. Using the local coordinates at each horizon we 
can now impose the Hamiltonian, momentum and Gauss-law constraints exactly as 


described in section |3.4| and obtain a set of generalised Stokes equations on each 
horizon. By solving these equations we can thus obtain currents on each horizon. 
Note that the precise source terms that appear in the Stokes equations on each 
horizon follows after imposing that E, ( are closed one-forms in the bulk and that 
they approach E^, at the AdS boundary. 

At the AdS boundary we can dehne the current fluxes through each d — 1 cycle Ca 
on Ed as in (3.49). Now consider any orientable d-dimensional manifold in the bulk 
with boundary Ca at the AdS boundary and a d — 1 cycle at the black hole horizon 
(which might be disconnected). Then using (5.3) and Stokes’s theorem, we deduce 
that the current fluxes are equal to the fluxes on the black hole horizon. In turn these 
fluxes can be obtained by solving the Stokes fluid equations on the black hole horizon, 
which only depend on the cohomology class of the sources at the horizon, which in 
turn only depend on the cohomology class of E^, on at the AdS boundary. Thus 
by expanding the fluxes E^X^ ^ basis of harmonic one-forms that are Poincare 
dual to the Ca we have a procedure for obtaining the DC conductivities. 

There are a number of interesting directions to pursue. In our analysis we assumed 
that the black holes have vanishing magnetic held. In [4] we will relax this condition 


and generalise the results obtained in 52,56 for Q-lattices. 

More generally, we have shown the DC conductivity of the boundary theory can 
be obtained by solving the generalised Stokes equations at the black hole horizon. In 
a narrow sense the huid equations are simply an auxiliary set of equations to solve 
this holographic problem. However, the innate physical character of the equations 
(with their novel viscous terms) suggest that their may be a deeper signihcance. As 
a hrst step, it would be interesting to determine whether the full time-dependent and 
non-linear generalised Navier-Stokes equations at the black hole horizon can also be 
used to obtain exact holographic information for the dual CFT. It is natural to expect 
that the time-dependent equations will be useful in extracting the small frequency 
behaviour of the AC conductivity. A related point is to use our results to develop a 
systematic hydrodynamic framework in the presence of holographic lattices. Finally it 
would also be very interesting to obtain some explicit lattice black hole solutions that 
depend on more than one of the spatial dimensions and analyse the fluid equations. 
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A Radial Hamiltonian formalism 


In this section we rewrite the equations of motion corresponding to the Lagrangian 
density: 


L = 


R - V(4,) 


mp2 


i (d-py 


(A.l) 


using a Hamiltonian decomposition with respect to the radial variable. We will 
follow the notation of [^, mutatis mutandis generalising to include the gauge-held 
and the scalar held. A useful reference is and we note that closely related work 


independently appeared recently in 59 


We consider a foliation by slices of constant r. We introduce the normal vector 
satisfying = 1. The D-dimensional metric induces a (D —l)-dimensional 
Lorentzian metric on the slices of constant r via hfj,^ = g^v — npi_ny. The lapse and shift 
vectors are given by = N{dr)f^ and — Nn^ where = {dr)^. In 

a local coordinates system we can write 


^ + jv“dr) {dx^ + N’^dr) , (A.2) 

where the shift vector has components = (0, A^“) and has components hrr = 
N^N’^'yab, Ka = labN^ and hab = Jab- Note also that = {N’^N’^'ybc'yabN'^)- 
We will decompose the gauge-held components via 

= ^ = -Nn^A^, (A.3) 

and hence = 6^ — N~^ $ In the local coordinates we have bj. = N°'Aa, ba = Aa 
and $ = —Ap -|- N°'Aa. 
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The radial Hamiltonian formulation can be obtained by first rewriting the La- 
grangian density as follows 

C =7Vv^ -V - 

- j , (A.4) 

where we have neglected total divergences. Here = \Cnhpy = pUy is the 
extrinsic curvature, K = g^’^Kpy and 

fdpby dybp , 

$ /d>\ 

+ —n VyUp — Dp > 

= fpyn^-^Dp^, (A.5) 

where the second expression utilises the fact that n^VpUy = —^DyN. This latter 
result follows from writing Uy = NVyT and using VpyV = XypV. We also recall 
that Dp is the Levi-Civita connection associated with the metric h and, for example, 
Dp<l> = hp’^Xy<l>. 

With hpy = Crhpy, bp = Crbp, 0 = Crcf) we can show that 

hpy = 2NKpy + DpNy + DyNp , 

K = ^^fp, + ^,iPpNP), 

= NnPfpp + NPfpp + XpibpNP) , 

0 = NnPdp(t) + NPdp(t). (A.6) 


The corresponding conjugate momenta are then given by 

SC 


ttP’' = = -V-h {Rp'' - K , 


flU 


TT^ 


Sh 
A r 

= ^ = y/DhhPPZXp , 


Sb 


= V^ZF- 


<fip 


n 


p ’ 




= —\/—hrFdy ((), 


(A. 7) 


where the second expressions for and which are not written in the canonical 
variables, are useful. 
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The Hamiltonian density, defined as "H = — £, can be written as a 

sum of constraints 


n = N H + N.H^ + ^C, 


(A.8) 


with 


H = - 


- 

^ D-2 
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- + ]^\f^hP^dp(t)da(t) , 

((_/,)-VV'^") + h^^f^pTiP 

- h'^^a^^/^Dp ((-h)-^/2 ^ h''^d^(j)TT^ , 

where tt = and we have ignored total divergence^ 

The equations of motion are given by 

hpy = -2N{-h)~^/‘^ (^Tipu-^T^hpu'^ +2D(^pN^), 

= -NV^ _ hd+i)j^j^„u ^ 

- ^ + 2Ar(-h)-^/2 

+ {D^D'^N - hP’' D^D^N) 

+ liV(-h)-'/2^-l (^TtPtT^ - hp-{hp^7lP7l-)^ 

+ liV^Z [hP^h-^hP^hpf^, - hp'^{hP^K<^fp^Us)^ 

- ^-N{-h)-^/^7:lhP'^ + (^PPh'^^dp4>d^4> - \{hP^dpct>d„ct>)hP'' 

+ v^T)^ ((-h)-^/2iVT'7r'^^) - 27r^('^T)^iV^) 


(A.9) 

(A.IO) 

(ATI) 


- VTTr 




d 


(A.12) 


as well as 


bp =Dp{N%,) - N (-h)-i/2 ^ j^u ^ 

{NZ h^PhP^fps) +2y/^D^ {{-h)-^/^ , 


(A.13) 


^In passing that we note that in the local coordinates ( A.2[ ) we have $ — N^ap = —A^ and one 
sees that using explicit coordinates one could uses Ar as a Lagrange multiplier instead of <&. 
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where we have dropped pieces proportional to the constraints in and We also 
have 


0 = - {-h) + N^dyct ), 

- N i V vr" + ^ /p. j . 

(A. 14) 


B Generalised Stokes equations from the constraints 


In this appendix we show how the Stokes equations (3.31 )-(3.33) arise from the con¬ 
straint equations (A.9)-(A.ll) in a radial decomposition. More precisely we will ex¬ 
amine the constraints for the perturbed metric, at linearised order, focussing on the 
leading terms of an expansion at the black hole horizon. In other words, we evaluate 
the constraints on a hypersurface of constant r and then take the limit r —)■ 0. 

We begin by noting that for the perturbed metric the unit normal vector has 
components 


„■ = g« n‘ = G-' (FU)-''^ 6gtr, 

= (B.l) 

Furthermore, the corresponding shift and lapse functions are given by 

= g"! Sgri, , 

The components of the extrinsic curvature take the form 

/(-« J- c-^u-^Ff-F'‘ t-a, (GU) + l^9r (GU) 6a„'\ 

^ 1 g-2j,-3/2p-./2 + a,(Gi7) AfJj , 

/f<. J_ 1 , 1 / 2 (Sg„-tGUQ'j +g-Ja, . 

K‘l = - u^l2p-i/2 + siSiSrSgt^ 

- U^FF-9^ s2‘9A95"'59™„ dr9iU ■ (B.3) 
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where here V is the covariant derivative compatible with the d-dimensional metric 
9dij • 

Expanding the extrinsic cnrvatnre close to the horizon we hnd 


^ — 


2 (47rT)V2 ^3/2 (^(0)3/2 


1 + 




GW 2 GW dvrT 


(0) -I 

+ —n*ailnG(°) 




2 (47rT)V2r-3/2 ^(0)3/2 
1 1 1 


(bp) 


(47rT)V2ri/2 G'Wi/2 






2 G(o)i/2r'/2 
and we also note that 


1 _ 1 


vW,, In G® 


2 GW 47rT 


,, ^ 1 (47rT)V2 


AttT 




Cit. 


(B.4) 


(B.5) 


2 G(0)l/2^1/2 

We now consider the following qnantity which appears in the momentnm constraint 


(A.IO) 


W, =D^ ((_/i)-V2^/.^) = + D,K 

= - ^ d^KxK^^ + d^K. 

Expanding at the horizon we hnd the following individnal components 

1 (47rT)i/2 X 


(B.6) 




2 ^{0)1/2 ^1/2 

111 


GWi/2 (47rT)V2 r-i 


17 ^ ( V(jni) - 27rT (j + - Vip 


where 


p = —2 tiT 


G(o) 


(rf+sS’) -■’ft-biSiVOnG'”’ 


(B.7) 


(B.8) 


Notice that after imposing the bonndary condition constraints (3.5) this dehnition 
of p is identical to the dehnition of pressnre given in (3.30). Another qnantity that 
enters the constraints is the momentnm of the scalar held. At leading order in r we 
have 




(B.9) 


We now tnrn to the gange held. From the second expression in (A.7) we have 


TT^ 


= y/^ZF^^nx = x/-gZF^^ = . 


(B.IO) 
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After expanding near the horizon we also hnd 


t 0, 

= fuJ' + ^ i^^w + E,) Z(°) af) (diW + E,) . (B.ll) 


We now now consider the constraint eqnations. Snbstitnting (B.IO) into the Ganss 
constraint (A.ll), = 0, we obtain the cnrrent continnity eqnation = 

0. When evalnated at the horizon this leads to = 0. We next consider the 


momentnm constraints (A.10) with lowered indices. Using the above resnlts we hnd 
that the t component gives = 0 while the i component gives the Stokes eqnation 

dPl. 


Finally, examining the Hamiltonian constraint (A.9), we hnd that the leading 
order expansion at the horizon for the terms involving the linearised pertnrbation 
implies Vjn* = 0, and hence gives no fnrther conditions. To see this we consider H in 


(A.9) as a snm of six terms and it is convenient to divide by y—h. Using (B.9) and 
(B.IO) we immediately see that the third and hfth terms vanish at linearised order. 
It tnrns ont that the leading order power of r that appears is r~^. We can show that 
the sixth term and, with a bit more ehort, the fonrth terms are of order r°. Next we 
consider the second term. The potential term is clearly of order r°. After examining 
the leading terms in the Christohel symbols we can also show that the Ricci scalar 
term is also of this order. Finally, we need to examine the hrst term. To do so it is 
convenient to note that nsing (B.4) we have 

1 


('-/i'l-iA n , _1_ I _ —V-v" 

1 ^ 2 (47rT)3/2G(0)3/2^3/2 


-h) 




2(47rT)V2G(o)i/2^i/2 

+ 4vrT^[o)(l - + 


1 




+ 




2G'{o) -rr ■ 4^2- • - ■ 

Continning to evalnate the hrst term we are eventnally led to the resnlt that the 
leading term in the Hamiltonian constraint can be written 

1 




-)■ 


2rG(o) 


V,;U. 


(B.13) 


C Holographic currents 


On-shell we have 

5S = j d^xV^ [4 5A;\ . 


(C.l) 
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where and are the radially independent, holographic stress tensor and current, 
respectively. After substituting the time dependent sources given in (3.1), (3.3) we 
hnd 


5S= + j\-tE{)\ 


(C.2) 


We thus see that —tEi is a source for the operator density G^f'^GWdj^ and —tQ is a 
source for the operator density G^!’^GWd^GG — /ij*). 

It is possible to show that the expectation values of these holographic tensor 
densities are given by 


G^/^TdiGG - ^If) = Q%- tG^/^VFdt'Xj • (C.3) 


Thus J*|oo and Q*|oo are the time-independent parts of the expectation values of the 


vector and tensor densities. To establish the first equation in (C.3) is straightforward. 


The second is a little more involved. Firstly, from the expression for Q* given in (3.20) 


and using (B.3) we can show that at linearised order we have 


Q* = EiGUf^G^, (^-2K^^ + - tGUQ^ - a,X . (C.4) 


Recall that if we write 


+ Xh^'^ + , 


(C.5) 


where X = 2K + /(0) -|- ... and Y corresponds to additional terms arising from 
the counterterms, then we obtain the stress tensor if we evaluate at the AdS 
boundary. Observing that at linearised order we have guh^^ = {5gtj — tGUQ) we 
can therefore write 

Q* = - YG - Ggtj - tGUQ)^ - atX . (C.6) 

We next want to take a limit as r —)■ cx). If the combination of R** and Y^^ that 
appear make sub-leading contributions, then we have 


Q% = hm [r^+^GY^G9~d (?' + FHQ) - 


(C.7) 


and we recover (C.3) after using that . We have explicitly checked for 

particular cases, eg H = 4 with Y^'' ~ that this does indeed occur. It would 

be interesting to hnd a universal argument that this is always true. 
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The time dependent piece on the right hand side of (C.3) is associated with the 


static susceptibility for the heat current two point function, as explained in appendix 


C of 36 . It can also be understood by noting that if we start with the background 
black hole geometries, with in particular t** = j* = 0, then the time independent 
linear perturbation that is generated by the coordinate transformation t —>■ t 
induces the transformation t** —>■ Promoting this perturbation to one that is 


linear in time leads to the time dependence as in (C.3). 


D Alternative derivation of the Stokes equations 


We discussed in section 3.3 how the generalised Stokes system of equations given in 


(3.31 )-(3.33) can also be obtained from the equations (3.9), (3.21) for J* and Q\ 


Indeed evaluating the hrst of the two equations in each of (3.9), (3.21) at the black 
hole horizon we immediately obtain = 0- These comprise two of the 


three Stokes equations, given in (3.31), (3.32). The third Stokes equation, given in 


(3.33), can be obtained from the second equation of (3.21). 


To obtain it we consider the pieces of Q* that are linear in r obtaining 




A 


- 4irTG<“'Cj 


( 0 ) ji 


(D.l) 


where we have dehned the matrix 


= 


— 


9 ( 0 ) 


4 

f2G(») 2G(»)y) 




(D.2) 


Notice that this matrix depends on next to leading order terms in the expansion at 


the black hole horizon. Equation (3.21) then implies that (D.l) should equal 


-d. 




G(0) 


— k I 


(D,3) 


To obtain an equation at the black hole horizon, we need to be able to express the 
matrix M in terms of leading order horizon data. After a long calculation, which we 
outline below, using the equations of motion for the background black hole we can 
show the key result 


= - 2 G (°) + 4^(0) 1/2 1/2 _ □^(0) 

+ ZG<»lal“>P« + , 


(D.4) 


and this leads to the hnal Stokes equation in (3.33). 
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We use the radial Hamiltonian presentation of the equations of motion for the 
background black hole solutions (2.3). The unit normal vector is n = dr- 

The lapse function is given by iV = the shift vector vanishes, = 0. 

The non-vanishing components of the extrinsic curvature are given by 


Ku = -^U^/^F-^/^dr{UG) , 

K,, = ]^U^I‘^F-^l‘^^rg^^, (D.5) 


and hence the non-vanishing components of the conjugate momentum are given by 


TT*^' = -if/ {Ggf''^ F-^/^g^’^g^^drOki + U^^^F-^/^g^^dr (UGgf^ , 

7i^ = {FG)-^/^g^/‘^Zdrat, 

vr^ = -GV2FVV/^0. (D.6) 


We also have 


IT = V’t"- = dr (UGg)'!^ 


(D,7) 


It is convenient to rewrite the equation of motion for given in (A.12) in the 
form 




- N _ (d+i)^ ^ y ^ {D^D'^N - D^D^N) 


+ -N^/^^h^^PF'^dpcIdd^ct) - hP^dp(t)d„(t)h^^'') , 


(D.8) 


where we used the fact that = 0 as well as the constraint i/ = 0 with FI as in 


dg). 

We now wish to plug in the background expansions (2.5) in the equations of 
motion (D.8). Taking the r —)■ 0 limit of the left hand side yields 

1 


_ 2N(-h)-'G A - i ^ ^ ^ (v^“’ . 

\/W) 


_ 2N{-h)-^/‘^ (TT*^ 7r\ - 5 TT TT*^' ) ^ ATAl + AnTgl^r , (D.9) 


( 0 ) 
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with as defined in equation (D.2). For the right hand side of (D.8) we find the 
following leading order behaviour as r —)■ 0: 


W ^ ('"‘’ 4 ) G'°’ - r,„,) 

+ . (D.IO) 


(ij) -> - G'°^g'4 ('^*'^(0) “ ‘‘*'^(0)SJJ) + ^(0) jJJ)) 

+ VG'(o) c/(o) 

+ (9;S)9g,9t'#>'°'a^'"> - 9 (o) 4^<”>9,9g,) ■ (D.ll) 

We now decompose the d + 1 dimensional Ricci tensor scalar via: 

\ V. ^ G-2 V.G V,G, 

= Q V. V.GWG^ GU, 

-^G'-2v,GWG. (D.12) 

We also have 

^ , (D.13) 

where V is the covariant derivative with respect to the d dimensional horizon metric. 
Putting these ingredients together we finally obtain the expression for given in 

ra. 


(rf+l)^ =(d)f>_ y 
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